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ABSTRACT 


Lindemann, Dieter. M.S.M.E., Purdue University. August 
1987. Experimental and theoretical studies on design cal- 
culations for latent heat storage. Major Professor: Prof. 
W. Leidenfrost, School of Mechanical Engineering. 

Research was conducted to obtain design data of latent 
heat-of-fusion thermal energy storage (LHTES) devices. The 
devices studied utilize a porous medium (PM) to increase 
. the effective conductivity of the storage system. Therefore 
the melting and freezing process is accelerated. Metal 


shavings and metal spheres were used as porous media and 


water as the phase change medium (PCM). 


A simple experimental method to measure the conduc- 
tivity of a water-metal mixture or ice-~metal mixture was 
applied to obtain conductivity data. These measurements 
can be considered fairly accurate because measured conduc- 
tivity of pure water and of pure ice agreed very well with 
those of accepted tables. The comparison of the measured 
conductivity for porous media with values obtained from 
several formulas from the literature showed that no formula 


predicts the conductivity with an acceptable accuracy. 


xviii 


Evaluation of defined quality values for LHTES' showed 
that metal shavings can enhance significantly the storage 


performance. 


As a practical application, calculation methods’ for 
LHTES systems were introduced and applied in a simulation 
program for residential space heating with energy storage. 
The simulation showed that metal shavings can help to save 
about 25% of the annual heating costs when compared with a 


LHTES system with pure water. 


CHAPTER 1 - INTRODUCTION 


In the long run, due to the limited amount of fossil 
fuels, we must depend on inexhaustible energy sources which 
require energy storage. Thermal energy storage is one pos- 
sibility which can be achieved in form of sensible heat 
storage or latent heat storage. Solid-liquid phase change 
is a form of latent heat storage. However its disadvantage 
is slow energy flow into and out of the storage system. 
This thesis investigates latent heat storage and how water 
saturated with metal shavings can be used as_the_ porous 
medium to enhance the energy flow. This is because the high 
conductivity of the metal can also increase the conduc- 


tivity of the water-metal mixture. 


“lel Exhaustible and Inexhaustible Energy Sources 


Exhaustible sources are coal, oil, and gas which are 
in limited supply. The inexhaustible sources are, on the 
other hand, continuously renewable. The earth” s biosphere 


receives energy naturally from three main sources: 


- intercepted solar radiation 


- thermal energy from the earth” s interior (produced 
largely by the decay of radioactive matter in the 


earth”s crust) 


- tidal energy 


Intercepted radiation can be found in application as: 
- solar energy 
- wind energy 
- hydropower 
- ocean thermal energy 
- biogas 


- wood 
Fige 1.1 shows how these sources can be used for Space 


heating. 


In early history men used almost exclusively inex- 
haustible energy. In this way, ships were propelled by 
wind, wagons by animals, and the houses were heated by 
wood, turf, or manure. Later, wind and water energy 
were applied for driving machines. In the beginning of the 
industrial revolution exhaustible energy usage began with 


the invention of the steam-engine powered by coal. 


Today we can already estimate the time when all 


exhaustible energy will be depleted. This is 


illustrated in [2] and shows that there will be no more 
U.S. production of petroleum liquids in about the year 
2060; therefore, we must now return to the use of inexhaus- 
tible energy sources applying modern standards’ of 


engineering knowledge. 


Exhaustible energy needs no special consideration for 
its storage because we find this energy already in a stored 
forme On the other hand, inexhaustible energy does_ need 
special considerations for its storage due to the fact 
that inexhaustible energy sources do not have a steady 
energy output. A time dependent energy resource can only be 
collected during a specific time frame. Our energy needs 
are widly varied and also time dependent but this depen- 
dance need not correspond directly to the collection time 


and therefore energy storage is required. 


1.2 Possibilities of Energy Storage 


Often energy has to be converted several times until 
it is finally used, and hence, there are many possibili- 


ties for energy storage. Let us consider the example of 


popcorn production from hydropower (Fig 1.2). In this 
example the energy can be stored in a reservoir, a bat- 
tery, or in thermal storage, or the product (here: pop- 


corn) can also be stored. It is our work to determine what 


the most economical storage possibility is. In this case, 


for someone who likes fresh popcorn, as opposed to someone 
who prefers stored popcorn. the choice is obvious when 
storage losses are considered. In situations where it is 
not so obvious what’ the “most economical storage possi- 
bility is, we must know more about the characteristic of a 
specific storage type. This thesis is concerned with the 


different possibilities for thermal energy storage. 


1.3 Thermal Energy Storage 


The major characteristics of a thermal energy storage 


are: 
- storage energy density 
- temperature range over which heat is added or removed 


- temperature difference associated with addition or 


removal of heat 
- temperature stratification in the storage unit 
- power requirements for addition or removal of heat 
- amount of losses from the storage system 


- cyclic durability of the storage 
The storage must be connected to other components of the 
energy systen. For this interaction between components, 


factors must be considered as follows: 


- period of time for which energy is stored or rate of 


energy losses 


- transportability of stored energy 
Heat may be stored by means of sensible heat or latent heat 
in a material. Sensible heat requires a temperature 
change; whereas latent heat requires a phase change or 
chemical reaction. The match between the energy source 
and the task, or load, is an important criterion in the 
selection of storage size and type. If large excursions 
in storage temperature are permitted, sensible-heat 
storage may be used. If in contrast delivery to load must 
be in a narrow temperature range, isothermal phase-change 
Storage using latent heat may be estautend. The type, 
size, and cost of a heat exchanger associated with the 
storage are closely related to storage thermal proper- 
ties. For example, a relatively small heat exchanger can 
be used if a high-temperature storage supplies a low- 
temperature demand. If the storage temperature is only 
slightly higher than the temperature demand, a large expen- 
sive heat exchanger will be required. If furthermore’ the 
energy is stored below the temperature of the demand 
system, a heat pump and hence even additional power is 


needed. 


1.4 Sensible Heat Storage 


~The best sensible heat storage materials have 
- high specific heat 

- high density 

- low cost 

- low toxicity 

- low inflammability 


- low chemical activity 


Sensible heat can be stored in liquids or in solids. 


leas] sensible Heat Storage in Liquids 


Water is the most common medium for storing sensible 
heat for the use with low- and medium-temperature systems. 
Water is gees and abundant and has a number of particu- 
larly desirable properties. The storage of water at tenmn- 
peratures above its normal boiling point requires expensive 
pressure vessels and is rarely cost-effective. Without a 
pressure vessel other liquids can be used for. thermal 
storage above 100 or Examples are Butanol, Octane, and 
Isopropanol with boiling points at 118 C, 126 C, and 


148 °c respectively. Unfortunately these liquids have 


lower density and specific heat with higher inflamma- 


bility and viscosity than water. 


1.4.2 Sensible Heat Storage in Solids 


For some applications sensible-heat storage in solids 
is advantageous. Hot air, for example, is often stored in 
Storage beds formed from solid particles like river rocks. 
The storage bed acts thereby as both a storage 
meet an and heat exchanger. Storage in beds has desirable 
thermal stratification capability arising from the 
poor particle-to-particle thermal contact in a bed of rock 
storage. It is rarely economical to use a single 
large mass of a solid material for sensible-heat' storage 
because a large mass has a very low surface-to-volume 
tatio, and hence the heat transfer from and to the storage 


is too small. 


Leo Latent—-Heat Storage 


Latent heat exchange can occur as solid-solid, solid- 
liquid, and liquid-gas phase change. Phase changes involv- 
ing gases should generally be avoided since the large 


volume of gas precludes compact storage. 


Solid-solid phase change usually involves’ hydration 


and dehydration. Na,SO, 10H,0 (Glauber”s salt) was the 


earliest phase-change storage material to be studied 


experimentally for house heating applications. At a tenm- 
perature below 32 ries the hydration point, the anhydrate 
becomes hydrate and crystalline with the evolution of heat. 
When the hydrate is heated above 32 °c the crystal dis- 
solves in the water and absorbs heat with a heat of 
fusion of 241 kJ/kg (compare with heat of fusion of water 
of 333 kJ/kg). The problem with salt hydrates is that 
the fusion temperature does not remain at the same _ ten- 
perature over several cycles. The cycle is not durable. 
After many cycles, the rehydration reaction requires 
progressively more subcooling to reach completion, thereby 


canceling the beneficial effects of isothermal phase 


change. 


Solid-liquid phase change materials should have sinmni- 
lar properties as those mentioned for sensible-heat storage 
materials above. In addition characteristics as follows 


should be present: 


high latent heat 

- an appropriate and well-defined melting point 

- high thermal conductivity 

- good tepeatability of the freeze-thaw temperature 


- low costs 


- low coefficient of thermal expansion 


- negligible phase change expansion 
If the material has not the last two mentioned characteris- 
tics, special care has to be taken for the design of the 
storage tank. Solid-liquid phase change materials generally 
have low conductivity especially in the liquid phase. 
As a result, energy flow into and out of storage is slow. 
This is also true for water. However, metal shavings in 
the water help to increase the conductivity of the storage 
medium because of the high conductivity of the metal. The 
storage system nearly keeps its energy density because the 
metal shavings have a high porosity which means that’ the 


volume of the phase change material nearly stays the same. 


1.6 Objectives and Strategy of the Study 


Objectives and strategy are structured into four parts 


and are: 


l. to find a method to determine the conductivity of a 
porous medium represented by metal shavings saturated 


with water or ice, respectivly. 


In this work a method is used for the first time 
where conductivities of porous media with water are 
measured by observing the water displacement due _ to 


freezing and melting. This measurement technique which 


aa 
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observes the integral effect of the phase _ change 
requires only to measure or to specify the temperature 
of the wall of the container in which the phase change 


occutsSe 


to define and find quality criteria (Chap. 3.4) for 
these porous media and conclude media with high 


performance. 


From basic equations quality criteria are deter- 
mined which show the capability of a medium to store 
energy and yield the energy density. With knowledge of 
the conductivity these quality values can be calcu- 
lated. From comparisons a porous media material of 


high performance can be found. 


to report about calculation methods for design of 
actual latent heat-of-fusion thermal energy storage 


(LHTES) systems. 


In practical applications mostly tube-in-shell 
heat exchangers are used. Equations will be given for 


a single tube as well as tube banks. 


to show the effects of a storage system with dif- 
ferent shavings in residential space heating applica- 


tion on the annual energy costs 


A computer simulation of a solar heated house 


11 


with a LHTES unit is used to show how different 


storage media can influence the annual heating costs. 


1.7 Background 


Investigations with heat pumps which use coils in the 
soil as the heat source have shown that during the winter 
ice is formed around the coils. This ice is formed in the 
soil saturated with water because the coolant temperature 
is below 0 °c. Low temperatures are necessary in order to 
achieve a high temperature gradient required to meet the 
energy demand. However, the frozen soil has to thaw during 
the summer, and problems arise when the summer is too short 
to thaw the soil. It has been proposed to use the ice 
built up by a heat pump in the winter for air conditioning 
in the summer. In this way the melting process would be 
completed in a shorter time and in addition surplus summer 
heat would be stored for heating in the winter. The melt- 
ing process can be enhanced furthermore by using materials 
of higher conductivity than soil, like water and >metal 


shavings as discussed above. 


Many studies have been undertaken to investigate the 
behavior of the solid-liquid phase change of water. 
Herrmann (1982) [3] and White (1984) [4] worked on “Melt- 
ing and Freezing of Water around a Horizontal, 


Isothermal Cylinder”, representing one type of a tube-~in- 
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shell heat exchanger. Weaver (1985) [5] studied “Solid- 
Liquid Phase Change Heat Transfer in Porous Media”. He 
worked with aluminum and glass beads in water. From his 
work there exists a cylindrical test cell and a program to 
calculate the freezing and melting process in the test 
cell. Both of these were used in the current’ research 
efforts. From the literature he reviewed he found five 
formulas which use different models to calculate the con- 
ductivity of a porous medium. Weaver found that these 
formulas give satisfying results for a porous medium con- 
sisting of water and glass but not for a porous medium 
with highly different thermal conductivities like aluminum 


and water. 
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Figure 1.1 energy sources for space heating. [1] 
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Example of popcorn production from _  hydro- 
power used to demonstrate the many existing 
storage possibilities for energy on the way 
from the source to its final use. 
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CHAPTER 2 - EXPERIMENTAL APPARATUS AND PROCEDURE 


The first objective of this thesis is to measure con- 
ductivities of porous media represented by metal shavings 
Saturated with water or ice, respectively. The conductivi- 
ties are measured by applying a melting and freezing pro- 
cess. This measurement method is unusual; nevertheless it 
is most similar to the melting and freezing in LHTES. The 
measurement techniques can be characterized like any other 


of the methods introduced in the next chapter. 


Ze Conductivity Measurement Methods 


Experimental study of the thermal conductivity was 
started in the eighteenth century. Since that time various 
methods have been suggested for different ranges of tem- 
perature and for various classes of materials having dif- 
ferent ranges of thermal conductivity values. Accordingly 
no method exists that would have optimum performance of 
conductivity measurements. Methods for the measurement of 


thermal conductivity fall into several categories: 


In the steady-state method of measurement, the test 


material is subjected to a steady temperature profile. The 
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thermal conductivity is determined directly by measuring 
the rate of heat flow per unit area and the temperature 


gradient. 


In the nonsteady-state method, the temperature distri- 


bution varies with time, and hence, the unsteady term in 
the heat equation, ~ —s is involved. With knowledge of 
density and specific heat the conductivity can be obtained 
from the diffusivity a. The nonsteady-state methods’ are 


further subdivided into periodic and transient heat flow 


methods. 


The heat flow has to be controlled in ae prescribed 
direction such that the actual boundary conditions agree 
with those assumed in the experiment. There are two possi- 
ble directions; and hence methods can further be subdivided 


into longitudinal and radial heat flow methods. 


There are two ways to obtain the rate of heat flow. 
In an absolute method, the rate of heat flow into a speci- 
men is directly determined, usually by measuring the 
electrical power input to a heater. In order to distin- 
guish between actual conducted heat and heat loss, the heat 
flow out of the specimen has to be measured. In a 
comparative method, the rate of heat flow is usually calcu- 
lated from the temperature gradient over a reference sample 
of known conductivity, which is placed in series with the 


specimen and in which hopefully the same heat flow occurs. 
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All methods are based on an experiment which can _ be 
modeled and mathematically described with high accuracy. 
An advantage is to have conditions in the experiment as 
they will occur in the application because additional 
errors resulting, for instance, from temperature dependence 


of the conductivity can be avoided. [6] 


Died Experimental Set-up 


The experimental set-up is used to measure conductivi- 
ties of porous media saturated with water respectively ice. 
It consists of the following three major components) intro- 
duced below: the test cell with heat exchanger, a constant 
temperature bath (CTB), and equipment used to observe the 
experiment. The integral effect of freezing and melting is 
observed by water displacement due to phase change expan-~ 


sion and contraction. 


The test cell, designed by Weaver, represents’ the 
model of a cylinder with an assumed adiabatic top surface 
and bottom and has uniform temperature on the circumfer- 
ence. The test cell is described in detail by Weaver in 
[5] and shown in Fig. 2.2.1. The capsule consists of a 
brass cylinder (158.8 mm long, 73.0 mm ID) which is capped 
on both ends by discs. They are pressed on the cylinder by 
five threaded rods and are made of acrylic plastic which 


has a low thermal conductivity (0.195 W/mK) which helps 
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reducing heat losses at the ends. Two joints are placed in 
the top cap to enable water flow in and out of the _ test 
section due to phase change expansion. Excess water leaves 
the enclosure through the joint in the center when the 
phase change material (PCM) freezes with inward interface 
motion. This water reenters through the joint near the wall 
when the PCM melts due to a wall temperature above freezing 


point. 


The counter flow heat exchanger consists of copper 
tubes wrapped and soldered around the outside of the brass 
cylinder which gives a uniform wall temperature. The test 
cell is insulated with Styrofoam fitted around the heat 
exchanger and the bottom cape The top cap is insulated by a 
piece of black insulation foam, which can easily be removed 
for observations of the phase change process.e- Most impor- 
tant is the insulation of the caps for the approximation of 
the theoretical model which assumes an adiabatic top cap 


and bottom. The additional insulation around the heat 


exchanger helps to avoid unnecessary heat losses. 


To supply the heat exchanger with cold or hot heat 
exchanger fluid, a constant temperature bath (CTB) 
Brinkmann-Lauda, RC 203 and another CTB Haake A82_ were 
used. The CTB circulates, with an internal pump, a water 
ethanol mixture consisting of 504 water and 504 ethanol 


through the heat exchanger. 


19 


Fig. 2.2.2 - 2.2.5 show the experimental set-up inclu- 
ding the measuring instruments. The main water displa- 
cement, due to phase change expansion and contraction, can 
be measured in the burette. During the process of freezing 
water moves from the opening at the center of the cell, 
through the capillary and into the burette. Similarly 
during melting water flows from the burette, through the 
capillary and into the test cell at the opening closest to 
the wall. The burette allows a reading error of 0.2 ccm. 
Smaller volume changes can be determined by observing the 
migration of an air bubble in the capillary. The air bubble 
is injected into the capillary with a syringe which can be 
stabbed through a septum. The septum is made from soft rub- 
ber and seals the puncture water-tightly after the syringe 
is removed. A septum is installed on either end of the 
capillary. The two septa on either end make sure that the 
air bubble can be moved to every location within the capil- 
lary by adding or removing distilled water with the 


sytinge. The capillary has an inner diameter D of about 


1.5 mm. The volume change V can be determined from the 
1 2 
traveled distance 1 of the air bubble by V = z * D4. 


The wall temperature can be measured with five copper- 
constantan thermocouples placed on the outside of the brass 


cylinder. 


In addition to the observations made with the burette 


and the capillary, direct visual observations are also 
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possible by looking through the transparent top cap of the 
test cell. When these visual observations indicate that 
the ice-water interface reached the center of the test 
cell, this fact is also indicated by a sudden change of the 
voltage of the thermocouples located at the center of the 
test cell. These observations are used to control and to 
verify the displacement measurements. Measurements of the 
temperature distribution in the test section are made by 
copper-constantan thermocouples arranged in four rakes with 
five thermocouples each. Radial and axial position of the 
rakes are shown in Fig. 2.2.1, position four and five. The 
thermocouples are connected with a digital thermometer 


which has a reading error of 0.l aol 


223 Conductivity Measurements with the Test Cell 


The test cell is used to perform experiments which can 
serve to measure conductivities for solid and liquid porous 
media. From PCM-freezing the conductivity of the solid 
porous medium (that means porous material with ice in the 
pores) can be determined. The wall -temperature and _ the 
freezing time have to be recorded. This method (compare 


with Chap. 2.1) is a 
—- nonsteady-state, transient 


- radial 
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- absolute 
method. The temperature distribution in the solid phase 
eoctes Sark time when heat is removed, and the solid-liquid 
interface moves inward. Liquid porous medium is still 
present in the center of the test section until the end of 
the freezing process. The heat flows radially outward 
through the solid porous medium, and the absolute amount of 
heat is directly dependent on the frozen mass and on the 
latent heat; hence measurements of power input are not 
necessary although the method is direct and not compara- 


tive. 


From PCM-melting the conductivity of the liquid porous 
medium (that means porous material with water in the pores) 
can be determined. The wall temperature and the melting 
time have to be recorded. The procedure is basically the 


Same as described for PCM-freezing. 


2-4 Standard Test Procedure 


In preparing for an experiment, the test cell was fil- 
led with the respective porous medium. Moreover distilled 
water had to be degasified by boiling for approximately 
fifteen minutes and then cooled to the fusion temperature 
while under vacuume This insured that the phase _ change 
process was not influenced by impurities from air bubbles 


in the test section. The water was then carefully siphoned 
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into the test cell with the porous medium already in it to 
avoid turbulent mixing which would enhance the absorption 


of air. 


The liquid PCM was then cooled down to a temperature 
of about QO °C using one CTB. In the mean time, the second 
constant temperature bath was cooled down to the freezing 
temperature selected for this experiment. The freezing 
experiment started when the pump of the CTB was switched on 
and the cold water-ethanol mixture started flowing through 
the heat exchanger. Freezing temperatures in the range of 
-25 "Cc to -5 °C were chosen. During the experiment, a 
representative wall temperature and the displacement’ were 
tead in intervals, depending on how fast freezing occurred. 
The monitoring of changes in the water level in the 
burette, due to phase change expansion, and in some cases, 
observation of an air bubble in the capillary, made it pos- 
sible to determine the time when all of the water in the 
PCM was frozen. In addition, temperature measurement in the 
test section and visual observations through the cap were 
made. From 25 available thermocouples, only ten tempera- 
tures were read with a digital thermometer. Those ten tenmn- 
peratures were two wall temperatures from representative 
Spots, five temperatures from the center thermocouple rake, 
and one temperature in the middle of each of the other 


three thermocouple rakes. 
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Similarly to the freezing experiment two melting 
experiments were carried out, one with horizontal and the 
other with vertical test cell position. The temperature of 
solid PCM was brought near to 0 °C using one CTB. [n the 
mean time, the second CTB was heated up to the melting tem- 
perature in question for this experiment. This CTB was 
then connected with the heat exchanger of the test cell to 
start the melting experiment. Melting temperatures in the 
range of 10 °C to 30 °C were chosen. Wall temperature and 


melting time were measured as described for freezing. 


2-5 Porous Media and Water/Ice Properties 


All relevant porous media properties are given in 
Table 2.5.1. For the experiment, beads and different 
shavings were used. In the “shape” column the diameter of 
the beads and the cross-section of the spirals is indi- 
cated. The value va is the ratio of volume and surface area 


of the porous medium. 
va =~ (2.5.1) 
a 


The aluminum beads are made from aluminum type 1100 with 
not more than 14 iron. The aluminum shavings were produced 
on a lathe from Al 6061 T4. Due to the different properties 
of aluminun, copper, steel, and brass it was not possible 
to produce the curls with identical geometry. Properties 


of water and ice are summerized in Table 2.5.2. and 2.5.3. 
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Figure 2.2.1 Schematic of test cell. [5] 
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Figure 2.243 Test cell resting on bottom insulation. 


Figure 2.2.4 


Insulated test cell with connection tubes 
constant temperature bath (CTB), digi 
thermometer, capillary, and burette. 
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Table 2.5.1 Porous media properties. 
11.D.# | metal Fig.# shape Va p Kn c.. ® I 
| l 
I mm mm kg/m W/mK JI/kak J 
ee ee 
i aluminum ~ beads diameter :3.2 0.533 2702 220 S80 O.424 | 
i aluminum = beads diameter:12. 2.0 2702 220 S80 0.532 | 
i 1 aluminum 2.5.1 wool 0.5x0.05 0.0231 2700 155 88o 0.939 | 
i 2 aluminum 2.5.2 chips 5.84x0.15 0.0743 2706 155 EEO o.903 | 
' 3 aluminum 2.5.3 curls 0.0508x5. 64 0.0252 2700 155 EEO 0.974 | 
! copper 2.5.4 curls 0.54.0 0.5 $933 419 3? 0.940 | 
| brass 2.9.9 797Cu, 3OAZn, curls 0.0484 BSS “106 375 O.956 | 
stainless tteel 2.5.6 AISI 304;curls 0.1x2 ©o.0464 7700 14.3 AE? 0.959 | 
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Table 2.5.2 


-1.0 


OWOMA NAW &F Wh we O 
ooooocoocqoco0o9$jo 


re) 
3 
kg/m 


998.156 
999.291 
999.447 
999.582 


999.696 


999.792 
999.868 
999.927 
999.968 
999.992 
1000.00 
999.991 
999.968 
999.929 
999.876 
999.806 
999.727 
999.632 
999.525 
999.404 
999.272 
999.127 
998.971 
998.803 
998.624 
998.435 
998.234 
997.077 
995.678 
994.061 
992.249 


Properties of water. 


of 
J/kg K 


4272.6 
4240.2 
4235.0 
4230.2 
4225.7 
4221.5 
4217.6 
4214.0 
4210.7 
4207.7 
4204.9 
4202.3 
4199.9 
O19 had 
4195.7 
4193.9 


4192.2 


4190.7 
4189.3 
4188.0 
4186.8 
4185.8 
4184.8 
4184.0 
4183.2 
4182.5 
4181.9 
4179.6 
4178.5 
4178.2 
4178.6 


9 
«x10 
mn /s 


124.24 
128.29 
129.03 
129.76 
130.47 
131.16 
131.83 
132.49 
133.13 
133.76 
134.37 
134.97 
135.56 
136.14 
136.70 
137.26 
137.80 
138.34 
138.87 
139.38 
139.90 
140.40 
140.90 
141.39 
141.87 
142.35 
142.82 
145.11 
147.28 
149.36 
151.32 


Table 2.5.3 


The 


satent heat. ot 
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p 
3 
kg/m 


920.300 
920.154 
920.006 
919.854 
919.699 
919.541 
919.380 
919.216 
919.047 
918.876 
918.700 
918.521 
918.337 
918.149 
917.957 
917.761 
917.560 
917.355 
917.145 
916.930 
916.710 
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Properties of ice. 


Cc 
J/kg K 


1961.3 
1969.1 
1976.9 
1984.7 
1992.5 
2000.3 
2008.1 
2015.9 
2023.7 
2031.5 
2039.3 
2047.1 
2054.9 
2062.7 
2070.5 
2078.3 
2086.1 
2093.9 
2101.7 
2109.5 
2T1L7 42 


a 
mn /s 


1.261E-06 
1.253E-06 
1.244E-06 
1.236E-06 
1.228E-06 
1.220E-06 
1.212E-06 
1.205E-06 
1.197E-06 
1.190E-06 
1.183E-06 
1.176E-06 
1.169E-06. 
1.162E-06 
1.155E-06 
1.1L49E-06 
1.143E-06 
.137E-06 
13 1E-06 
.125E-06 
1.119E-06 
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Figure 2.53.1 
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Aluminum chips. 
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Aluminum curls. 


Figure 2.5.3 
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Figure 2.5.5 





Brass. 
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CHAPTER 3 - ANALYSIS 


The resistance method is one of many which can be used 
to calculate freezing and melting. This method is applied 
to calculate the freezing and melting time in the test 
cell. The conductivity can then be calculated from the 


predicted and measured time for freezing and melting. 


i es | Calculation Methods 


All methods used for measurements of thermal conduc- 
tivity are based on an experiment which can be modeled and 
mathematically described. The first step for this descrip- 
tion is to calculate the freezing and melting process which 
is done by calculating the fusion front position as a func- 


tion of time. 


An exact analytical closed form solution can not _ be 
obtained for the geometry of freezing or melting in a 
cylinder because nonlinearities introduced by the boundary 
conditions preclude exact solutions ina closed domain. 
However, in cases like this two methods commonly used yield 
a solution: the integral method and the finite difference 


method. 
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"The application of integral methods to the _ solution 
of problems involving the solidification of liquids ini- 
tially at fusion temperature" is shown in [7]. Here for- 
mulas for a cylinder are given. The "Karman-Pohlhausen 
method" [7] assumes a simple one-parameter temperature 
profile and yields a simple formula for the freezing time 
as a function of the interface position. Unfortunately 
this method appears to fail when applied to the circular 
cylinder. The "Tami method" [7] assumes ‘a two-parameter 
quadratic profile which makes the solution much more 
involved requiring numerical integration. This method still 


has only an accuracy of about 10 &. 


The finite-difference method is applied by Weaver in 
his program for the cylinder. A finite-difference method 
is accurate to the extent allowed by the assumptions on 
which the method is based. Therefore the finite-difference 
method is used here to show the accuracy of the simple 


resistance method. 


3.2 The Resistance Method 


The resistance method [8] analyzes freezing and melt- 
ing on the basis of a variable resistance analog to account 
for the motion of the interface and the varying thickness 


of the solid or liquid phase, respectively. 
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Fig. 3.2 The resistance method uses a variable resistance 


analog. 


3.2.1 Calculation of Freezing and Melting Time 


Fig. 3.2.1 shows a schematic diagram of the test cell. 
The cylinder is uniformly filled with porous medium and 
liquid initially at uniform temperature of 0 "Gis The fusion 
front starts moving radially inward when a uniform wall 
temperature a less than fusion temperature Ty for freezing 


or T > T. for melting is imposed on the cylinder inside 
Ww 


surface. 


The assumptions made in the calculation are as_ fol- 


lows: 
O 
1. Initial temperature is 0 C. 


2. The porous medium is isotropic (properties are not 


dependent on directions). 
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3. The porous medium is homogeneous and properties are 


constant in space. 


Ls Conduction is the dominant mode of heat transfer in 


both the solid and the liquid. 


Ds There is no bulk fluid motion due to the difference in 
density between the solid and liquid phase and natural 


convection is absent. 
6. The porosity is constant. 


7. The local temperatures of the phase in the voids and 
the porous medium are the same ( iee., there is no 
heat transfer between the porous medium and void phase 


constituent). 


8. The fusion front movement is quasi steady and the 


temperature profile changes only slightly with time. 


9. The sensible heat is negligible compared to the latent 


heat ( h. >> c [T,-T ] ). 
f P f ow 


Freezing will be developed first. With the above assump- 
tions we have a temperature distribution in the solid 
porous medium as in a hollow cylinder with a thermal resis- 


tance due to conduction 
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_ _in (R/r) 


Re cond  2ani1k Coney 
sm 


where 1 is the length of the cylinder and et the conduc- 
tivity of the solid phase of the porous medium. When the 


contact resistance at the interface formed by the wall and 


the ice layer is neglected, the total resistance is 


1 
= Lec 
an l i (3 
en + aE 
R 
t,cond t,conv 
but with R =o in a solid we get 


t,conv 


R 
t t,cond 
Hence the heat rate flowing through the solid is 


toot 
q=- > | (3.253) 


Re cond 


where the minus sign indicates heat release. The heat’ rate 
leaving the cylinder can also be calculated from the latent 
heat multiplied by the rate of mass which is undergoing the 


phase change 


Sip ee 
a” TE dt 
dr 
q os he Piso 27TYr 1 dt ® (3.2.4) 


Note that h, is defined for the density of ice OP; Since 


f , O° 
O 

the interface is at 0 C, the density has to be determined 

at this temperature. Substituting Eq. (3.2.1) into Eq. 


(3.2.3) and realizing that Eq. (3.2.3) is equal to Eq. 
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(3.2.4) yields with separation of variables and integration 


h,. 9 ® 
i (T,-T ) dt = - tA. for 1n(R/r) dr 
sm 


The wall temperature is often not constant in real applica- 


tion. For instance in the experiments for conductivity 
measurements (see Chap. 4) a transient wall temperature was 


observed. From the integration of the left hand side we get 


T_-T dt = T. dt-o- T dt 
f (t-1) J 
= T t - T t 
w,a 
* (T, 7 Tae 


where T is defined as the average wall temperature 
w,a 


l 
oe = t f TT dt C3 6225) 


After the integration is performed the time t required for 


freezing to a radius r can be found. 


Do 
>| 


( = r° InR - r~ [ = Inr - = ] ) + C 


The integration constant is obtained from the initial con- 


dition 


r(t=0) =R 


which yields the integration constant 
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with which finally 


P, 
f i,0 \ 4 Py2 2-2 
= (OO 7 rr eoLe 
t= oper Fim r/R) + G(R -r™)) (3.2.6) 
sm f£ wa 
or with 
r/R =x C3624-7) 
2 
QO, ® R 
£ Pico 1 2 1 2 
wer a is eLe 
t (TAT y 7 x Inx 61 x )) (3.2.8) 
sm f 'w,a 
The total freezing time t. = t(r=0) for freezing of the 


whole cylinder volume is obtained for the limit of t as r 


goes to zeroe With the law of 1° Hospital 


2 
limr In(€r/R) = 0 


r>Q 
he Pp, 0 ® R° 
t = limt = ee (3.2.9) 
: r>Q : Kom g : ) 


The derivation for melting is similar except that the 
heat of fusion has to be transported through the liquid 
porous medium. In contrast to the freezing case we have for 
the melting case convection as well as conduction heat 


transfer. The total resistance is again 


1 
Ra Ge een (3.2.2) 
Ma, Be eee 


R R 
t,cond t,conv 


Unfortunately R is unknown for a porous medium in an 


t,conv 


enclosure of arbitrary shape. Therefore we define 
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| in(R/r) 
. = 


= (3.2.10) 
t 27 1 K oFe 
Where k ff is an effective conductivity including both the 
e 
influence of conduction and convection. However K fF is 


also not known yet but will be obtained from experiments. 
In contrast to the conductivity K om for the solid porous 
medium, the effective conductivity Kore depends on the spe- 
cial geometry. We now follow the same procedure as done 
already for freezing. Heat has to be supplied to the system 
and hence we get a change of sign in the equation 
corresponding to Eq. (3.2.3) where heat was released. In 
the equation which corresponds to Eq. (3.2.4), po is still 
the density of ice at 0 “iO Following the same steps 


further we obtain 


+e dateey & -(R°-r")) (3.2.11) 


ct 
il 
n 
“of 
a | 
ar 
_w 
om 
ho 


h. 0, ® R 
t =< at = 5 x* Inx + ¢(1-x7)) 012) 
eff-ow,a f 
hee. Os 6 R 
= f *~ 1,0 
A ea CT | (362643 


3.2.2 Check of the Resistance Method with the Finite- 


Difference Program 


Weaver” s program [5] for freezing in porous media in 


a cylinder was used to check the resistance method. The 


46 


assumptions on which this program is based are assumptions 
2 through 7 from Chapter 3.2.1. The program solves the heat 
equation for the solid and liquid phase, respectively, 
using boundary conditions at the wall, the interface, and 
the center line together with the initial condition for the 
temperature oof the porous medium (PM). The problem is that 
the location of the interface is not known a priori; there- 
fore the heat equation has to be solved for small time 
steps repeatedly as the interface moves inward until the 


whole volume is frozen. 


For the check of the resistance method, freezing of 
water (porosity: = 1.0) at an initial temperature 


O 
T, = Q C in the test cell was calculated. The results are 


shown in Table 3.2.2. 


Table 3.2.2 Comparison of freezing time for resistance 
method and finite-difference method for 
water initially at 0 sare 





Freezing Time 


CF Resistance t. Finite-Difference 
Method Method ih 


43.6 min | 46.4 min 6.0 


ee 
174.2 min 175.4 min 0.7 
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The error of the resistance method seems to be smaller 
or at least comparable to the error of the Tami method [7] 
(about 10 %) which requires many more calculations. The 
error of the resistance method is small if the wall tem- 


perature is close to 0 °C because assumptions 8 and 9 match 
much better with reality. The total freezing time calcu- 
lated from the resistance method is always less than the 
freezing time from the finite difference method because the 


former neglects the fact that the ice must be cooled down. 


3426.3 Conductivity Obtained from Average Wall Temperature 
and Freezing Time 


The instantaneous wall temperature is known from. the 
experiments at certain intervals and the freezing, respec- 
tively melting time is recorded. The average wall tempera- 
ture ae is calculated from Eq. (3.2.5) which can practi- 
cally be done with time weighted averaging, trapezoidal 
rule, or Simpson” s 1/3 rule. The conductivities can now 
be calculated from Eq. (3.2.9) and Eq. (3.2.13). For freez- 
ing 
Pe ee Cs a (3.2.14) 


and for melting 


f i 
i 
Keft 4 t_(T =1.) eee) 
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363 Approximating the Conductivity for the Liquid Porous 











Medium 
In contrast to the effective conductivity K fF is the 
true conductivity k not influenced by convection. This 


lm 


conductivity can be obtained from a first order approach 
based on two melting experiments: melting in a vertical 


cylinder and melting in a horizontal cylinder. 


It is possible to calculate the conductivity ie from 
k and a correction factor F. 
eff 
= CPs ee 
Kan K vee F ( ) 


We define the time ratio 





c= = (3636529 
E 
t,v 
In this equation is t the total time for melting in the 


t,h 
horizontal cylinder and Peg the total time for melting in 
the vertical cylinder. The convection influence is of 
course different for these two geometries and therefore 
also the melting times will differ. The correction factor F 
depends on the amount of convection and this can be seen in 


the time ratio ec. We assume F = F(c) to be linear. Two 


points of this function are already known. 


First, for the time ratio c of water F must be chosen 


so that F corrects K fe to Kin = SH0” the conductivity of 
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water which is known from tables. Hence 


k 
H,0 
ECCS 6? & 
2 eff, 0 
Z 
Second, for c = 1 no convection is present (if there 


would be convection there would also be a difference 


between horizontal and vertical total melting time) and 


thus 


F(c=1.0) = 1.0 


Now the function F(c) can easily be drawn. 





Fig. 3.3 The correction factor F as a function of the 
time ratio c. 


1.0 if c > 1.0 
F(c) = A ctB it Cy oo €< c <1.0 C3343) 
2 
a 
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where 


k 
HO 
a= — (3.3.4) 
Seg 
HO 
A = — (25335) 
HO 
l as 
H,0 


3.4 Quality Criteria for Porous Media 


There are several requirements for energy storage in a 
porous medium. These are previously mentioned in Chap. 1.4 
and 1.5. Some requirements are influenced especially by 
the liquid in the PM; others are influenced by the PM 
itself. This thesis gives its attention to the PM and con- 
Siders only water as the material undergoing phase change 
in the eer medium. Three new dimensions will be defined 


to quantify the storage behavior of the PM. 


The heat storage power P shows the capability of a 
medium to store heat ina freezing or a melting process. 
High energy transfer Q in a short time gives a favorable 
high power P. As _ reference value we take the energy 
transfer necessary for total freezing of a volume and_ the 


related total freezing or melting time. 


Q 


p= — (3.4.1) 


t 


a 


Q. = (Ie p (1-8) a ee °. 6] T +h ®) V (3.4.2) 


t ice ice m f Piio 


with volume V, porosity $6, and mean temperature 
l 
T = > (T + T_) (3.4.3) 


For the cylinder we can take the results from the resis- 
tance method for the total freezing or melting time ti 
Eqe (3.2.9) or Eq. (3.2.13), and substitute them together 
with Eq. (3.4.2) into Eq. (3.4.1). 

k 


f Pio 
4n1 ABS(T._-T ) (3.4.4) 
f Ww 


®] T +h 


pS ({c 0 (1-8) ae ice m f Piso 


p 


ice h ® 


The first line in this equation is a first quality value, 


the performance factor 


k 


(3.4.5) 
f Pi 0 


pf=([c p (l-®)tc, o, @)]) T +h d ) 


ce ice m f Pio h © 


The performance factor pf is independent of geometry which 
is described with the geometry factor (here for example 
4 ~ 1 for the cylinder) and independent of the driving 
potential ABS(T, - TO. ABS is the absolute value function. 
The conductivity k is the solid porous media conductivity 
an for freezing and the effective conductivity K fe for 
melting. For phase change in a sphere the geometry factor 
would be 8 nm R where R is the radius of the sphere, and 


2 A/x would be the geometry factor for the plane wall, 


where A is the area of the wall and x the fusion front dis- 
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tance from the wall. The performance factor has units of 


conductivity W/mK. 


The energy density Q°”” is a second quality value. 


= [cp (1-8) +c, 6) T +h ® (3.4.6) 


icePice m £ Pio 
Since we are looking for a number showing the overall 


suitability of a porous medium, we multiply pf and Q°”” to 


get a number which we call the storage factors 
s = pf oe (3.4.7) 


For the evaluation of the experimental data a. = -10 C 


10 °C for melting. 


was taken for freezing and T 


3.5 Calculated Conductivities of Porous Media 


For a system with a single phase occupying the void 
and a single porous medium, different models can be found 


to calculate the overall thermal conductivity. 


The series model assumes that the porous’ medium 
behaves like a set of alternate layers of fluid and solid 
perpendicular to the mean heat flow. This assumption gives 


the thermal conductivity 


ie. sop eaten (3.5.1) 


where the subscript lm, m, and 1 refer to the effective 
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liquid property influenced by the porous media, the porous 
medium and the liquid phase occupying the void, respec 
tively. Similarly, the corresponding equation for a solid 
phase occupying the void is obtained by merely replacing 
the subscript 1m and 1 with sm and s, respectively, which 
refer to the effective solid property influenced by the 


porous medium and the solid phase occupying the void. 


The parallel model also concerns a set of layers, but 
the mean heat flux is parallel to the direction of the 


layers. 


Kin = kT ay kK, ® C3654:2) 


The Lichteneker geometrical model is listed in [9] 


i: 
2-8-(1-a)1/3 4 = ((1-3)1/3-(1-8)) 
Nc i, ee ee 
Lm 1 K 
1 - (1-@)1/3 4 = (1-3)!/3 


m 


Also listed in [ll] is the empirical model which can be 


obtained from dimensional analysis. 


® 1-@ 
k =k x! 


ln 1 mn (3.5.4) 


Veinberg [10] proposed a model which claims to be univer- 
sally applicable for a medium with randomly distributed 


spherical inclusions. The nonlinear implicit equation is 


ag DD 


m = 
1 
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In [11] is the Zehner and Schluender model set out. 


k 
lm 2 1-@ 
ete { l- 1-6 ] + 1=1G 
1 
& . + = 
(1 ADE in(=4) _ Sth 2 “= (3.5.6) 
(1-26) | 
where 
1-6_10/9 
G=H[ 3 ] 
and 
KI 
nN = i 
m 


The geometry factor H is 1.25 for spheres and 1.4 for  bro- 
ken. particles. This formula is not defined for all combi- 


nations of ky? k and ® because for 


sie 


1-6,10/9 


k, 7 1625 ky [937] 

the denominator 1-AG becomes zero. For this special case a 
formula is given in [11]. In order to illustrate the first 
five formulas and to show the different behavior of the 
functions, these formulas are plotted for a water metal 
System and an ice metal system, respectively. The plots are 
three-dimensional plots showing the conductivity K and the 
porosity ® in the horizontal plane and the result, the 


effective PM conductivity, on the z-axis vertical. The 


related Figures are Fig. 3.5.1 through Fig. 3.5.10. 
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The parallel and the series model are the limits of a 
possible effective PM conductivity. The other models men- 
tioned above have different shapes lying embedded within 


these limits. 


Comparison of experimental data with the Zehner and 
Schluender formula shows that the formula overpredicts for 
shavings by about 50 % and underpredicts for the beads by 
about 60 4 In the following text it is not further 
accounted for the Zehner and Schluender formula. The vali- 
dity of the other formulas is checked in Chapter 4 and 


Appendix A. 
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T » interface 





solid in porous media 


Figure 3.5.1 Schematic of the test cell. Freezing with 
inward interface motion. [5] 
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Figure 3.5.2 





Thermal conductivity of water and porous 
medium calculated by using the series model. 
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Figure 3.5.3 Thermal conductivity of ice and porous 
medium calculated by using the series model. 
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Figure 3.5.5 
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Figure 3.5.6 Thermal conductivity of water and porous 
medium calculated by using the Lichteneker 


geometrical model. 
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medium calculated by using the Lichteneker 


geometrical model. 
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Figure 3.5.8 Thermal conductivity of water and porous 
medium calculated by using the empirical 
model. 
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Figure 3.5.10 
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CHAPTER 4 - EXPERIMENTS 


Previous to measurements of the conductivity of porous 
media, experiments were performed only with distilled water 
in order to check, the apparatus, the method, and the accu- 


racy of the measurements. 


4.0.1 Experimental Data 


Thirty-five experiments were performed for the thesis. 
The data of these experiments were evaluated by a computer 
program. The output of this program is given in the Appen- 
dix A. Each experiment covers two pages in the Appendix. 
The first page describes the porous medium, melting or 
freezing conditions, total time for phase change, measured 
conductivities, conductivities for comparison, especially 
those obtained from five conductivity models, and values of 
the three defined quality criteria. On the second page a 
plot of the phase change as a function of time is shown. 
The solid line shows the observed phase change _ process. 
Every star on the line marks one measured value. The 


dashed and dotted lines show the phase change process’ cal- 


culated with the resistance method and average wall 
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temperature from Eq. 3.2.8, Eq. 3.2.12. The percentage of 
phase change is calculated from the radius of the interface 


by 
ae ; 
percentage of phase change = (l-x ) 100 2% 


The dashed line uses the true conductivities without con- 
vection influence. Only for melting a dotted line is plot- 
ted which shows the phase change as it would be with a con- 
stant effective conductivity for the whole experiment. The 
dotted line for melting and the dashed line for freezing 
always lie above the solid line because the actual wall 
temperature in the first part of the experiment is closer 
to 0 "6 than the average wall temperature. This means that 
the actual driving potential in the beginning is smaller 
than assumed, because using the average wall temperature 
throughout the whole experiment yields less phase snenee: 
For comparison, values of interest are taken from Appendix 


A and are given in tables incorporated in this text. 


4.0.2 Determination of Beginning and End of Phase Change 








For accurate measurements, the exact freezing and 
melting time has to be determined which means that the 
beginning and the end of the phase change has to be ascer- 
tained exactly. This was done by observing the movement of 
the air bubble in the capillary. After the flow of the 


coolant through the heat exchanger was initiated for a 
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freezing experiment, the air bubble first moved quickly 
towards the test cell because the test cell shrank when the 
wall was suddenly cooled. Thereafter, the same movement 
continued when the water cooled down. When the ice forma- 
tion started, the air bubble jumped towards the  burette 
where it left the systems. However, when freezing experi- 
ments having an average wall temperature of about -10 °C or 
less were performed, it was observed that ice was already 
formed a few seconds after the coolant flow was initiated. 
This means that for low temperatures, observations of the 
air bubble were not necessary. The same behavior was 
observed for melting experiments - the melting started 
nearly without any delay after the flow of the CTB-fluid 
through the heat exchanger was’ started. The end of the 
phase change, that is the time when everything in the test 
section is frozen or thawed, was more difficult to detect. 
The percent volume change in a cylindrical test cell is a 
logarithmic function of time which means that the volume 
change rate gets smaller and smaller with time, making it 
more difficult to determine the end of the phase change. Of 
course here also observations of the air bubble achieved a 
higher accuracy compared to observation of the water level 
in the burette. Observations with air bubble and _ burette 
were checked by temperature measurements in the test sec- 


tion. When a thermocouple makes contact with the fusion 


front, its temperature changes within a few seconds. The 
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thermocouples on the center rake were therefore used _ to 


check the phase change end. 


4.1 Measurements with Distilled Water 


o<e lel Freezing 


Freezing with temperatures of the coolant of less than 
-10 °¢ instantly gives an ice layer on the inner cylinder 
wall. In the vertical test cell position this ice layer 
grows almost exactly cylindrically inward. For horizontal 
freezing however, a water volume shaped like ae funnel 
remained unfrozen below the top cap for. about four minutes 
so that the end of the freezing process was much less shar- 
ply defined. Thus the vertical position of the test cell 


was chosen for later freezing experiments. 


Dendritic ice forms when water is supercooled consid- 
erably below 0O °C before nucleation occurs. It consists of 
crystals which form thin plates dispensed in the water. The 
ice growth occurs only in the supercooled region and ends 


O 
at the 0 C isotherm. The quantity of the dendritic ice 


that is formed can be related to the amount of subcooling 


m, | O (T,-T.) 
LCe- 2 

m h 
H, 0 f 


where TH is the nucleation temperature. 
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Subcooling and formation of dendritic ice was observed 
in the freezing experiment #3 with an average wall tempera- 
0 ee eer O 

ture of -4.1 C. The water initially at -0.1 C cooled 
e ® O 

down within 43 minutes to a uniform temperature of ~-4.0 C. 

Initiated by a hit on the test cell, nucleation occurred. 

Dendritic ice was instantaneously formed, shooting from the 

wall right into the center of the test cell since all _ of 
0 e 

the water was at atemperature below 0 C. The water in 

O 

between the dendrites was again at 0 C after the dendritic 

ice was formed. The burette showed a displacement due to 

nucleation of 3 cem. Thereafter normal ice formation 

started at the cylinder wall, incorporating the ice den- 


drites, until the whole test section was frozen. 


Subcooling and dendritic ice formation should _ be 
avoided for LHTES systems because it precludes’ heat 
transfer at a constant temperature. Furthermore, it 
requires a longer freezing time because water which con- 
ducts heat before nucleation has a lower conductivity than 


ice. 


4.1.2 Melting 


For melting, none of the two positions of the test 
cell yield exact uniform radial interface movement because 
of convection effects, meaning that the applied phase 


change model describes the melting process not fully exact. 
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In the vertical position, the warm water moves to the _ top 
so that the ice at the top of the test cell melts faster. 
Accordingly, at the end of the experiment the rest of the 
ice is approximately shaped like a parabolic cone placed in 
the middle of the bottom of the test cell. For horizontal 
melting, it could be observed that an ice rod remained in 
the center of the test cell diminishing in diameter. The 
melting did not occur concentricly but aroynd a line below 
the center line of the cylinder, due to the fact that’ the 
rising warm water causes a higher melting rate on top of 
the rod. At the end of the experiment the center of the ice 
rod was approximately 10mm away from the center line of the 
cylinder. Nevertheless, when the remaining piece of ice 
finally lost its contact with the thermocouple rake, it 
floated up and was quickly melted due to direct contact 
with the warm cylinder wall. The ice rod in these melting 
experiments showed grooves cut into the ice-water interface 
in planes perpendicular to the axis of the cylinder. The 
grooves were deepest under the ice rod and became shallower 
towards the top. In experiment #6 with an average wall ten- 
perature of 25.7 "GC. five grooves were present. The 
grooves show the presence of free convection which is, in 
an enclosure, based on convection cells, where one convec~ 


tion cell belongs to each groove. 


4.1.3 Freezing and Melting with Initial Temperatures Dif- 


ferent from Fusion Temperature 


Three experiments were performed in which the initial 
temperature was different from 0 aoe Table 4.1.3.1 compares 
freezing and melting time of these experiments with the 
equivalent ones which had an initial temperature of about 
o °c. Comparison of experiment #2 and #9 shows that freez- 
ing of water with an initial temperature of 19 °C took 94% 
longer. Also the finite-difference program was used to ver- 


ify the measured freezing times of experiment #2 and #9. 


The finite-difference program calculated a total freezing 
time th of 67.7 min for experiment #2 and 76.4 min for 
experiment #9, hence a percentage difference of 13% which 
is a slightly higher difference as calculated from the 


measurements with 9. 


Table 4.1.3.1 Comparison of freezing and melting time for 
experiments with initial temperatures equal 
and unequal 0 aor 





reezing m: melting 
v: vertical h: horizontal 
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For melting, respective values of 11% and 20% were 
obtained. These percent values of the time difference were 
in the range of the values obtained from Senge the 
ratio of sensible and latent heat. For freezing in a 


°C) are calculated 


cylinder, values ee t.(T,)/t, (7, =0 
using Weaver” s finite-difference program [5]. The calcu- 
lated time ratio for the cylinder is independent of size. 
These values listed in Table 4.1.3.2 can be used to correct 


Freezing times obtained from methods which assume the ini- 


tial temperature to be the fusion temperature. 


Table 4.1.3.2 Ratio of total times for freezing in a 


cylinder - freezing from an initial tempera- 
ture greater than 0 7G divided by freezing 
from a temperature T, = 0 C. 
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4.1.4 Conductivity Measurements for Ice 


Table 4.1.4.1 Comparison of measured conductivities of ice 
with conductivities from tables. 


E T z 
*P W,a ice table 
if mC W/mK W/mK 
1 vertical -18.6 Z2«16 Dis 21 263 
2 vertical -12.0 2.18 2.20 0.9 
3 vertical -4.,1 1.76 2.18 19.3 
4 Leo 


Experiments like the freezing experiment #3 where sub- 
cooling occurred are not suitable for conductivity measure- 
ments. Those experiments have a longer freezing time com- 
pared to a similar experiment without subcooling. Hence the 
measured conductivity from experiment #3 is too low. The 
other conductivities measured in experiments with average 


O 
wall temperatures of less than -10 C show only an error of 
about 2 % when they are compared with conductivities from 


tables evaluated at the mean temperature from Eq. 3.4.3. 


4.1.5 Conductivity Measurements of Water 


Table 4.1.4.2 shows the measured effective conduc-— 
tivity of water calculated from Eq. 3.2.15 which is not a 


property of the medium but rather depends on the geometry 
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1 


Table 4.1.4.2 Effective conductivity of water in the test 
cell averaged over total melting time. 












5 vertical 26.8 
6 horizontal Za 22-64 
7 vertical 9.6 1.00 
8 horizontal 9.5 


of the vessel. The effective conductivity is not even con- 
stant for one experiment because the enclosure gets bigger 
when melting proceeds, and so the effective conductivity 
increases. Thus a low effective conductivity is present in 
the beginning of the melting process and a relatively high 
conductivity is present in the end. Since there is more 
convection heat transfer in the horizontal cylinder than in 
the vertical cylinder, the effective conductivity is higher 
in the horizontal cylinder. The experiments show also that 
higher convection activity is achieved with higher wall 
temperature from which also a higher effective conductivity 
results. The measured true conductivity (without convection 
influence) for water is the exact value because the calcu- 
lations from Chapter 3.3 for the true conductivity uses 
just the known conductivity of water to calibrate the equa- 


tion. 


Nevertheless, there is still a possibility to calcu- 


late the true conductivity without the convection influ- 
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ence. For a Rayleigh number less than 1000, there is no 
convection in a cylindrical vertical or horizontal enclo- 
sure. From knowledge of the Rayleigh number, the maximal 
gap width of the molten region can be calculated for which 
the test cell has only conduction . The result shows that a 
water layer of 3.4 mm can be molten in the test cell 
without having convection. An effective conductivity calcu- 
lated with Eq. 3.2.12 from the phase change right in the 
beginning in the gap hence yields the true conductivity 


k, . The true conductivity k Fe is a property of the mediun. 


lm 1 
Table 4.1.5 Comparison of measured conductivities of 
water - without convection influence - with 
conductivities from tables. 
| Exp T k Ra=1000) k Az, 
W,a H,0 table 
i °C W/mK W/mK y 
6 horizontal 2561 0.493 0.604 -18.4 
7 vertical 9.6 0.615 0.574 7.1 


Table 4.1.5 compares measured true conductivities of water 
with conductivities from tables. The conductivity is cal- 
culated for a gap width equivalent to that of Ra=1000. The 
percent error for water here is higher than the error for 
ice because here only values from the very beginning could 
be used. In the beginning the change of the wall tempera- 
ture is high (see Chapter 4.1.6), and the reading of the 


wall temperature was limited to time increments of 30 
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seconds. On the other hand, the gap width of 3.4 mm_ was 
already reached after 2 to 4 minutes which means that there 
was not enough data present to calculate the average wall 
temperature (Eq. 3.2.5) with a higher accuracy and there- 


fore, the accuracy of the true conductivity Kon is limited. 


4.1.6 Wall Temperature 


The plot in Fig. 4.1.6 gives an example for variation 
of the wall temperature with time for experiment #1 - vert- 
ical freezing. When freezing starts, only a thin ice layer 
separates the fusion front from the wall, which means only 
small resistance for the energy transfer and high load, and 


O 
hence, a wall temperature close to 0 C. When the ice layer 
grows, the load decreases and the wall temperature con- 


verges towards the CTB temperature. 
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Figure 4.1.6 
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temperature of the constant temperature bath 
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4.2 Freezing and Melting in Porous Media 


Freezing and melting were observed in the test cell 
for different porous media in order to measure their con- 
ductivity and to evaluate their performance for heat 
transfer enhancement, and usefulness for energy storage. 
Experiments with different kinds of aluminum porous’ media 
were performed to get information about the influence of 
the shape of the PM on the phase change. Also, the 
influence of the material was studied by performing exper- 
iments with copper, brass, and stainless steel in addition 


to the experiments with aluminum. 


From temperature measurements with the thermo couples, 
the shape of the fusion front was observed. The fusion 
Front for freezing in the test cell is cylindrical 
throughout the whole experiment. For melting, the same 
deviation from the cylindrical shape was observed as obser- 


ved already for pure water. 


Normally a clearly formed interface was observed. 
However, for the 12 mm aluminum beads and _ the copper 
spirals, the metal pieces reaching out from the interface 
act as fins. For example, during freezing, the main inter- 
face moves gradually inward while ice also builds up around 
the beads or spirals, acting as fins. The ice that builds 
up around these fins is later incorporated into the main 


ice mass when the main interface has moved further on. For 
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melting, similarly, it was observed that the ice thawed not 
only around the cylinder wall but also around the spheres 
and spirals which are in contact with the wall of the test 
cell. It seems that for this type of phase change, there 
has to be a high conductivity of the PM present, as_ well 
as a high ratio of volume to surface area, va, as in big 
spheres or heavy spirals. To which extent this type of 
phase change influences an overall conductivity measurement 
and calculation for a geometrically straight surface could 


not be determined. 


Dendritic ice formation is not as likely in porous 
media as it is in water because the porous media yield more 


possibilities for nucleation than pure water does. 


es | Freezing with Different Shapes of the Porous 





Medium 


The porosity of optimal packed spheres in an infinite 
container is 0.26. Nevertheless, the beads with 12 mm 
diameter from experiment # 24 and the beads from experiment 
# 12 with 3.2 mm diameter have a higher porosity because 
the wall effect prevents the optimal arrangement for 


highest density of the beads. 


The conductivities of the aluminum porous media 
increases with decreasing porosity. This is predicted by 


all five formulas from Chapter 3.5 which are used _ (for 
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Table 4.2.1 Values obtained from freezing in aluminum 
porous media. 


Exp ® k K Qu" pf S 
sm emp 
iA “107° “107? 
: 3 4, 
it W/mK i, I/m W/mK JW/m K| 
beads Des 8.5: +6 § % . 6 ~ 28 
3.2mm 


LE RED eT TD ED cI enn SMS 





calculations and it is also evident since aluminum has a 


much higher conductivity than water. Amazing however is the 
big difference in conductivities of wool and curls although 
their porosity and volume to surface area ratio is almost 


the same. 


For all experiments the empirical formula (Eq. 3.5.4) 
predicted the conductivity best. (Only for the beads with 
3.2mm diameter the Veinberg model calculated the conduc- 
tivity more accurate.) The absolute value of the error is 
between 15 and 68%. The highest energy density has, of 
course, the medium with the highest porosity because the 
major fraction of stored heat is stored in the form of 
latent heat which is achieved only in water and not in the 
aluminum. The performance factor, pf, follows the trend of 
the conductivity and is of same magnitude. The storage fac~- 


tor, s=Q°°~’pf, obtained by the beads is the highest. This 





| 
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happens due to the fact that even a small increase of the 
metal volume causes a relatively big increase of the con- 
ductivity; whereas the energy density varies linearly with 


the porosity. 


Ge 2s Melting with different Shapes of the Porous Medium 


Table 4.2.2 Values obtained from melting in aluminum 
porous media. 





i W/mK: W/mK W/mK W/mK Iw/m'K 

13/14 beads 7.81 7.81 6.81 6.81 1.3/7 
3.2mm 

25/26 beads 6.89 6.89 7243 7.43 1.63 
12.mm | 

19/20 chips 3.04 1.70 4.01 1.94 1.45 

22/23 curls 3.07 2.42 3.60 Zeal. . 1.39 

16/17 wool 22-02 0.68 2.95 0.68 1.16 


v: vertical 
h: horizontal 


When melting data is discussed, it has to be recalled 
that conduction and convection are of importance. The true 


conductivities Kan for water and aluminum are all lower 


than the conductivities ae for ice and aluminum because 
the conductivity of water is only about one fourth of the 
ice conductivity; however, the convection influence is also 
of importance. Convection in addition to conduction, which 


are both included in the term k increases this value in 


eff’ 


comparison to the true conductivity, k For the beads, 


ln 


however, the melting times in the vertical and horizontal 
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test cell position are the same, hence it follows that con- 
vection has no effect, and effective and true conductivi-~ 
ties are the same. The effective conductivities follow the 
assumed pattern - the conductivity is high when the poro- 
sity is low. The convection influence can be seen in the 
plots in the appendix when the solid line is compared with 
the dotted line. Convection is present when the _ phase 
change rate is low in the beginning but high in the end 
(nearly no decline of the slope in the solid curve) com- 
pared to the dotted line. The ratio of K EE divided by Kin 
shows the convection influence, where a high ratio signi- 
fies high convection activity. This ratio is for curls, 
chips, wool, and pure water 1.3, 1.8, 2.9, and 3.0 respec-— 
tively and 1.3, 2.1, 4.5, and 4.5 respectively for the hor- 
izontal position. The wool allows almost as much convection 
as pure water. Although the chips have a smaller porosity 
than the curls, convection in the chips is of greater 
importance. In the horizontal test cell position more con- 
vection is present. The effective conductivity of ice and 
aluminum is about two times the conductivity of water and 
aluminum. The storage factor for freezing is of the same 
magnitude as the storage factor for melting, only for the 
beads is it in the melting case considerably lower. In all 
but one of he. ween experiments, the best results were 
obtained again by the empirical model for sonduceteees cal- 


culation. 
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4.2.3 Freezing of Water in Different Porous Media 


EEC GREECE ee 


Material 
Table 4.2.3 Values obtained from freezing of water in 
copper, aluminum, brass, and steel porous 
media. 
Exp k K om aa Q pf Ss 
A% hor i | 
: 3 ; 4 
if W/mK W/mK Ze J/m W/mK JW/m K 
2/7 copper 409.8 5.79 -48 0.307 6.19 1.90 
21 alumi- 155.0 3.93 -38 0.317 4.2 1.32 
nun 
33 steel 14.3 2.44 = 0.313 2.60 0.81 


In addition to copper, brass, and steel, the aluminum 
curls from experiment # 21 are included in this table 


because they have a comparable shape and porosity. 


The conductivity of the saturated porous media 
increases with increasing conductivity of the porous media. 
The porosity was about the same in all of these exper- 
iments. The empirical formula yielded the best results. The 
energy density is nearly the same for all of these exper- 
iments since the porosity does not vary a lot. The copper 
porous medium achieved the best results in respect to con- 


ductivity and storage factor. 
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Ae 25h Melting of Ice in Different Porous Media Material 


mp «EEC eo ee 


Table 4.2.4 Values obtained from melting of ice in cop- 
per, aluminum, brass, and steel porous 
media. 








E K kK h: K hi 
oe eff Lm Wage ee i 
r W/mK W/mK W/mK W/mK  JW/m°K 
28/29 copper 2.83 1.40 Sek 1.51 1.40 
22/23. “aluni-= 3.07 2.42 3.60 2.71 1.39 
num 
31/32 brass 2.19 1.14 2.80 1.22 1.07 
34/35 steel 1.89 0.88 2.69 1.04 1.03 


vi vertical 
h: horizontal 

The conductivities increase with increasing conduc- 
tivities of the porous metal. Higher than expected, 
however, are the conductivities of aluminum. Copper 
achieved the best storage factor. The ratio of effective 
Over true conductivity is about the same for all materials 
with the exception of aluminum which shows a lower ratio. 
The best predictions of the conductivity mere obtained by 
the empirical model. The values of the empirical model were 


again too low. 


Gs 2D Improvement of the Empirical Formula for Conduc- 





tivity Calculation 


As shown above the empirical formula is the best 
method for calculating true conductivities. Although this 


formula reasonably follows trends in variations of porosity 
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and conductivity of the PM, the obtained values are almost 
always too low. Therefore a correction factor of 1.35 is 
proposed to represent similar conditions as those investi- 


gated in this thesis. In this way it is obtained 


® 1-6 
= 1.35 k, _ ) Cae 2 654-1) 
sm ice m 
_ ® (1-0) 
Kon 135 “H,0 Kin (4290-3462) 
for 0.7 ¢< 6 ¢ 0.98 


In this way only an error of + 35 & has still to be expec- 
ted. For the given porosity range, a lower limit for the 
conductivity can be obtained from the original empirical 
formula Eq. (3.5.4) and an upper limit, although mostly far 
off, from the parallel model Eq. (3.5.2). Eqe (4.2.5.1) 
and Eq. (4.2.5.2) should only be used for the given poro- 
sity range because due to the correction factor the known 
conductivities for @ = 0 and @ = l1 are not’ regained 
correctly anymore. It could not be found out how the ratio 
volume to surface area, the length of the spirals, and 
phase change with the fin effect affects the true conduc- 
tivity. Nevertheless it seems that these influences are 
small compared to variations in the metal conductivity and 


especially in the change of porosity. 
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CHAPTER 5 - FREEZING AROUND A SINGLE TUBE AND IN 
TUBE-IN-SHELL HEAT EXCHANGERS 


LHTES units normally consist of tubes in a storage 
tank through which the coolant flows. Such an arrangement 
is called tube-in-shell heat exchanger. In this chapter 
first calculations for freezing around a single tube in an 
infinite shell are presented which are then extended _ to 


calculations of tube-in-shell heat exchangers for LHTES. 


The calculations from Chapter 5.1.2, 5.2.1, and 5.2.2 
are also valid for melting with turbulent free convection 
Flow. For this case k has to be replaced by _ k where 

sm eff 
K ret has to be obtained then from experiments with a scale 
model. This is possible because the convection coefficient 


in turbulent buoyancy driven flow is independent of the 


characteristic length. 


del Freezing Around One Single Tube in an Infinite Shell 


Pe ew | Exact Solution from Ozisik 


For freezing around one single tube in an infinite 


shell an exact solution is given by Ozisik [12]. This solu- 
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tion accounts for initial temperatures of the liquid dif- 
ferent from the fusion temperature. The solution given by 
Ozisik accounts only for conduction. Since initial tempera- 
tures of the liquid which are different from the fusion 
temperature are allowed, convection is present in the 
liquid phase. The convection in the liquid phase can be 
accounted for approximately when in Eq. (5.1.1.2) the true 
conductivity k is substituted by the effective conduc- 


lm 


tivity k The calculation is only valid for a line heat 


eff. 
sink which could be approximated by a heat exchanger tube 


with a very small diameter. The position of the fusion 


front has to be calculated from the transcendental equation 


Dole Lie 1s 
2 k (T -T_) ane / 
Q -n’ Im i “sm “1m _ 42 
7S + ————_~_—__———- e = "a QO, h 
Tt Ei(-h-a@ Ja. ) sm “i,0 f 
sn lm 
where Q is the heat sink strength in W/m and 
k 
det =a C54.) 5 lix-Z.) 
m pc)i., 
= + l- Oe aol we 
(pc), Ppic, (1-®)p_c ( ) 


Similar formulas are obtained for a and (pc) if the 

sm sm 
subscript 1 is replaced by se The similarity parameter A 
is 


i. (5.1.1.4) 


2(a_ etl? 
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-Ei(-x) is the exponential-integral-function which is 


defined for x > O as 





-Ei(-x) = f = dt (S215) 
l 


values for the exponential-integral-function are given in 


Table 5.1.1. 


5.1.2 Resistance Method 


The resistance method can also be used to calculate 
Outward freezing around a cylindrical tube. The same 
assumptions have to be made as for inward freezing in a 
cylinder (Chapter 3.2.1). 

ee Ot 2 


1 R,1 2 252 
t = k T —T ) [5 LE In=+7(r —-R ) ] (5.1.2) 
sm ff w 


3.2 Freezing in a Tube-in-Shell Heat Exchanger Consisting 


of a Tube Bank 


3-2.1 NTU Method 


Shamsundar et. al. [13] calculated freezing in a 
tube-in-shell heat exchanger with the NTU method. In these 
calculations is accounted for temperature variation of the 
coolant in the axial direction but axial heat conduction is 


assumed to be negligible. Further assumptions are: 
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- PCM initially at freezing temperature 

- cocurrent flow 

- tube geometry as shown in Fig. 5.2.2 

- sensible heat negligible 

The method accounts for forced convection in the tubes) and 
requires therefore the Biot number which can be calculated 


from 


Bi = CSe2e151) 





where D is the tube diameter, k is the measured conduc- 
sm 

tivity of the solid porous medium and h is the convection 

coefficient calculated from correlations for internal 


forced convection. It is a dimensionless time defined by 


k t 
O h. o D2 0 f in 
Pio f 
T is the inlet temperature of the coolant. Shamsundar et 


in 


al. state that this definition and calculation method 
accounts for a time dependent inlet temperature. Another 
necessary value is obtained from 
NTU = n DL h/(m Cc), (502% 13) 
The fraction F of PCM that is frozen in a_ radial plane 
O 
element at the inlet to the storage can be calculated by 


iteration from 


F 


_ O 
ec. = J f Bi + OCF) ] (5426144) 
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with the auxiliary function OCF) given in Fig. 5.2.2 and 


J = ~ - — for aligned tubes 


or 


for staggered tubes. 


co |e 


The effectiveness € is defined as 


ee (5.2.1.5) 
; T .-T. 
f in 
or 
Ee = aes: eee (54.2416) 
m c (1. -T, ) 


where ma is the specific heat of the coolant and € can be 


obtained from iteration of 


NTU = -ln(l-e) + Bi [Q(F )-Q[Cl-e)F 1] CS 20147) 


where QCF ) is another auxiliary function given by 
Fige 5.2.2. The energy transfer rate Q can now be obtained 
from Eq. 5.2.1.6. The frozen PCM fraction F can finally be 


obtained from 
F = (Foe + Bi [¢(F )-¢(1-e)F\])/NTU (5.2.1.8) 


This method also applies to the counter-current flow 
heat exchanger over the first phase of heat recovery during 
which the tubes of the heat exchanger have no. mutual 
interaction. When the tube pitch is not close to the tube 


diameter, this first phase forms the bulk of tthe freezing 
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time. During the remaining period, the counter current heat 
exchanger will be a better performer than the cocurrent 
onee Using the data pertaining to the latter will hence 
enable a conservative estimate. For other ratios of pitch 
and tube diameter as shown in Fig. 5.2.2 two integrals have 
to be solved. The input to these integrations is obtained 
as the result of two-dimensional finite-difference calcula- 
tions in a representative plane perpendicular to the _ tube 
axes. The details of these calculations can be found in 


[13]. 


5.2.2 Resistance Method for Tube-in-Shell Heat Exchanger 


As long as the fusion fronts of the different tubes do 
not make contact with one another the resistance method 
derived for a single tube is valid. When the fusion fronts 
then finally meet, a kind of inward freezing can be assumed 
to calculate the last part of the phase _ change. In this 
chapter the calculation for a staggered tube bank of the 


design indicated in Fig. 5.2.1 is shown. 


After the fusion fronts make contact, further inward 
freezing of the remaining unfrozen shape is estimated to 
take place like freezing in a circular cylinder. This vir- 


tual cylinder has the radius 


= 2 






OY 
Ge 
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Fig. 5.2.1 Freezing around staggered tubes 
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th ut radius of heat exchang tub F 
f the freezin g time fron given f n £ tion 
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it must be distinguished if a certain amount of solid phase 
(and hence frozen fraction F) will form before or after the 


fusion fronts meet. The frozen fraction F at the time when 


the fusion fronts meet is 








l 2 2 
_ 4 L - R 
er = ae ; ; (oS, 2e262) 
- R 
27 
When F < F é the interface radius r can be calculated 
mee 
from the given frozen fraction F with 
1/2 
3i/2 2 = 


the freezing time t can now be calculated from Eq. 5.1.2. 


When F > F : the radius of the interface in the virtual 
mee 


cylinder is 


1/2 
op Oe a 3 (5.2.2.4) 
Vv An 2 ae 

and 

t = t(F ) + [t€r )-t(r .)] (5.2.2.5) 

meet Vv v,i 

Be is the radius of a cylinder with the same volume as 

b ) 


the volume of the liquid phase in a triangular unit when 


the fusion fronts meet. Therefore r |. = r (F ) from Eq. 
v,i v meet 


(5.2.2.4) which yields 





-5) = 0.11328 L 
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From the equation for outward freezing 


h. 0 D 
= Ly, _ ke re! Rete 2Rol op be 22 
6 t(>) = or 7 [oh In7—+;7 (Gh Red) (5<262.6) 
sm f w 
from the equation for inward freezing 
h. Pp. ® r 
f "i,0 1:2 are ree eee 
a ne sal foes Panik al 
t(r_) RUT TS [sr in(s )+F (RY ri] (5.2.2.7) 
sm f w V 


and also from the equation for inward freezing 


h. 0 ® r 
£.. 24,40 1 2 Vote age 2... 2 
BF is, ee pe en = 
ECE a? C=) oa R ie tae © NeeeeeS) 
The Eq. 5.2.2.2 / 5.2.2.3 / 5.2.2.4 are derived in 


Appendix B. 


Table 5.1.1 Values of the exponential-integral-function. 








x E,tx) x E,(x) x E,(x) x E(x) 
E(x)or — En — x) Function 

| 0.00 x 0.25  1.0442826 0.50 0.5597736 1.60  0.0863083 - 
0.01 4.0379296 0.26 0.0138887 0.55 0.503364! 1.65  0.0802476 

[12] 0.02 3.3347078 0.27 0.9849331 : 
0.03 2.9591187 0.28 0.9573083 0.60 04543795 1.70  0.0746546 
0.04 2.6812637 0.29 0.9309182 0.65 04115170 1.75 —0.0694887 
0.05 24678985 0.30 0.9056767 0.70 0.3737688 1.80 0.064713! 
0.06 2.2953069 0.31 0.8815057 0.75 03403408 1.85  0.0602950 
0.07 2.1508382 0.32 0.8583352 
0.08 2.0269410 0.33 0.83610!12 080 0.3105966 1.90  0.0562044 
0.09 1.9187448 0.34 0.8147456 O.X5 0.2840193 1.95  0.0524144 
0.10 1.8229240 0.35 0.79421S4 0.90 02601839 20  4.89005( — 2) 
0.11 1.7371067 0.36 = -0.7744622,--0.95.-0.2387375) 2.1 4.26143 
0.12 1.6595418 0.37 0.7554414 
0.13 1.58%8993 0.38 0.7371121 1.00) -0.2193839 2.2 3.71911 
O14 14241457 0.39 0.7194367. 1.05 0.018728 2.3 -3.25023 
0.15 14644617 - 040 0.7023801 1.10 01859909 24 284403 
0.16 1.4091867 0.41 0.6859103 1.15 0.171558 2.6) 2.18502 
0.17 1.3577806 0.42 0.6699973 
0.18 1.309796! 0.43 06546134 1.20 O.1584084 2.8 1.68853 
0.19 1.2648584 0.44 0.639732K 1.25 463134 3.0 1.30484 
0.20 1.2226505 045 = 0.6253313 1.30 O13S84S10 3S) &9TOTS( = 3) 
0.21 1.1829020 $0.46 @.6113865 1.35 O2SS1OK §=—4.0 3.779358 
0.22 1.1453801 0.47 -0.5978774 1.40 0.1162193 4.5 2.07340 
0.23 11098831 O48 058847843 148 01077774 5.0 1.14830 
0.24 1.0762354 0.5720KKK 1.50 


ee we eee +: 


0.49 


0.1000196 


wv. 


0 


(The figures in parenthesis indicate the power of 10 by which the numbers to the 
left, and those below in the same column, ure to be multiplied). 
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(94) h 





storage 


Auxiliary function for phase change 
calculations. [13] 


Figure 5.2.2 
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CHAPTER 6 - APPLICATION OF LATENT HEAT STORAGE 


Finally, the usefulness of a storage system depends 
not only on quality criteria like the ones defined and dis- 
cussed above but also on the way a storage tank is incor- 
porated in a system and how the parts of this system act 


together. 


In order to find out more about these interactions, a 


System was modeled and analyzed by a computer program. 


6.1 The Model 


Fig. 6.1 shows the schematic of the model. The storage 
tank serves as a seasonal storage in space heating applica- 


tion. Cooling was not considered. 


The storage tank is a cylinder with two heat exchanger 
tubes, se: from the collector and one from the heat pump. 
The tubes are situated on the centerline of the cylinder. 
Calculations for the storage tank are performed with the 
resistance method. In order to avoid freezing of the 
Storage tank from the outside during the winter, the tank 


is buried and placed at a location with a distance z to the 
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surface, where zis the largest freezing depth of the soil 
during the year. In this way the storage tank even 
receives heat from the soil; an effect which is also known 


from heat pumps with ground coils. 


The heat losses of the house are calculated with the 
degree-day-method [14]. A loss coefficient-area 
(UA) = 463 W/K was chosen. It is assumed that the house 
has no capacity at all to store heat not even for short 


periods. For the heat pump an efficiency of 0.3 was chosen. 


Collector calculations were performed as shown in [14] 
with collector data from Fig. 7.7.1 of [14]. Meteorological 
data were taken also from [14] Appendix G for a site in 


Madison, Wisconsin. 


For the control system of the house it was decided 
that water from the collector is passed through the radia- 
tors of the house whenever the temperature of the water 
leaving the collector is higher than a minimum temperature 
of 35 °¢ and heating of the house is necessary. If the tem- 
perature of the water from the collector is less than the 
minimum temperature, the inlet to the radiators at A (Fig. 
6.1) is closed and the heat is used for energy storage (by 
melting the ice in the storage tank) as long as the tenm- 
perature of the fluid is above 0 °C. Whenever the collector 
power is not enough for space heating, the use of the heat 


pump is prefered for additional heat supply rather than 
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using external energy from a power plant or a conventional 


furnace. 


The program simulates the heating system during one 


year. It starts calculating in April and assumes that the 


Storage tank is filled by a frozen fraction F = 0.99. 


Approximately the same number was obtained at the end of 


the year. Iteration is necessary to calculate the elements 
of the system because the calculation result from one ele- 


ment has an impact also on performance calculation of other 


elements. 


Parameters used for the simulation were as follows: 


Site Madison,Wisconsin ; 

latitude 43.1 

climate 4294 degree-days 

per year 

House loss coefficient-area 463 W/K 

(UA) 

mine required 35 °C 

collector fluid tempera- 

ture for heating 

O 

Collector slope 63.1 

azimuth angle 0.0 0 

surface area 40 m 

flow rate 0.133 kg/s 
Storage tank length 50 m 

diameter 0.5 m 

depth z 2.0 m 
Heat pump efficiency 0.3 


i, 


ES LES OE 
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6.2 Results 


The energy necessary for heating the house over. the 
period of one year is Q = rs Pe ae kWh. Table 6.2 shows the 
influence of different porous media on the performance of 
the system. This influence is documented by the ratio of 
auxiliary energy and total energy necessary to heat’ the 
house. Auxiliary energy Oiace is needed to drive the heat 


pump or to heat the house directly. Table 6.2 shows further 


values of the quality criteria for horizontal melting. 


Table 6.2 Results from simulation program compared 
with quality values for horizontal melting. 





PM —Q 1002 Qe" = pf Ss 
° -9 ® -9 | 
10 10 
“| 
yA J/m W/mK JW/m K 
without storage 85 | 
only water 74 0.348 3.0 1.05 
CK ¢¢ 2.64 W/mK) 
Saturated porous | 
metal | 
steel 71 0.335 3.07 1.03 
Al wool 69 0.345 3235 1.16 
beads 3.2 66 0.161 8.47 Ls37 
brass 62 0.334 342 1.07 
Al curls 61 0.340 4.1 1.39 
beads 12. 61 0.196 8.95 1.76 
Al chips 57 0.317 4.59 1.45 
copper 55 0.329 4.25 1.4 
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The simulation showed that the storage tank was too. small 
for the beads. Already in the early summer the PCM in the 
tank was thawed and no more heat storage was possible. In 
case of the shavings, the PCM was thawed in September or 
November. When pure water was used, a frozen fraction of 
about 0.1 remained. The problem was to thaw the frozen mass 
from the previous winter with a reasonable collector size. 
Table 6.2 shows that quality criteria can give hints for 
Storage design. However, the simple multiplication of 
energy density, 0O°°”°, and performance factor, pf, yields 
results which might put too much emphasis on the perfor- 


Mance factor. 





Power Plant 


Latent Heat Storage 


Figure 6.1 Schematic of the modeled house for  simula- 
tion of latent heat storage in real applica- 
tion. 


COT 


104 


CHAPTER 7 - CONCLUSIONS AND RECOMMENDATIONS — 


jem | Conclusions 


No formula was found which predicted the conductivity 
correctly, but a simple apparatus was described to measure 
the conductivity. 

Shavings used as porous media enhance the performance of a 


storage tank because they 


O 
keep the nucleation temperature TD close to 0O Cand 
avoid subcooling even at freezing temperatures only 


slightly below fusion temperature 


- increase the apparent conductivity of the liquid and 


solid phase 


- increase the performance factor pf and _ therefore 


increase the added or removed heat rate 


- reduce the energy density only very little due to the 


high porosity. 


From the experiments it can be concluded that the optimum 


porosity © for a LHTES with water as PCM should be between 
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0.6 and 0.9, where a value of about 0.6 is appropriate for 
high heat rates that have to be added or removed. A value 
of about 0.9 should be chosen when energy density is more 
important than heat cranes ets A porous medium with a high 
conductivity yields the best results. However, the conduc- 
tivity of the water-metal or ice-water mixture does not 
increase proportionally with an increase of the metal con- 
ductivity. Hence it can be economical to use a cheaper 
metal and to increase the heat transfer rate by other 


meanSe 


When the conductivity is known from any experiment 
freezing in a tube-in-shell heat exchanger can be calcu- 
lated from the NTU-method. For melting, on the other hand, 
the effective conductivity has to be Pac. een experiments 
with a scale model; and this procedure is only valid for 


turbulent flow in the melt region. 
7-2 Recommendations 


Further research should address three problems: 


1. Find an accurate formula to predict the conductivity 


of a mixture consisting of a PCM and a porous medium. 


2. Find a calculation method for the melting process’ of 


ice filling the voids of a porous medium. 





36 
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Perform experiments with tube-in-shell heat exchangers 
and determine the accuracy of the NTU-method and the 


resistance method. 


The procedure for 2. would be as follows: 


Calculate heat transfer due to free convection in an 


enclosure filled only with the PCM (6=1.0). 


Account for the obstructions of the porous’) medium. 


(Introduce the permeability of the porous medium). 


Account for the enlargement of the liquid region when 
melting progresses and its influence on free convec~- 


tion. 


Account for the irregular shape of the fusion front 
which is formed due to convection when melting 


proceeds. 
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Appendix Ae Experimental Data 


Explanation to Appendix A are given in Chapter 4.0.1. 
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Experiment # a 
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Experiment # 4 
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Experiment # ? 
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Experiment # 3 
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Experiment # 9 
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Experiment # 10 


metal > no metal only DISTILLED WATER 
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Experiment # 14 
metal > ALUMINUM 
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Experiment # 1& 


metal >  ALUINUM 

porosity . phi = 0 FOS 
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density 1 oT ts cere F700. 0 kg/m#te 
specific heat: Ch = 390.0 J/kg 
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faperiment # if 
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# 22. ALUMINUM_CURLS 
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no eb .00 25.00 320.00 25.00 46.00 45.00 5.00 
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Cxperiment # 24 


diameter: Leimm 


ce. OOOO mim 


metal > SLUM IMUM 
porosity - pha = 
description : HEAIMS, 
¥Gl. “area MEL. 
densaty 7 can 


specific heat: cm = 


FREEZING i4: ea VEN PIC 


2702.0 kasmtte! 
880.0 J/kQgkth 


AL cylinder 


temperature GY beth oo... ee. Tb = ~-14.90 
ambient temperature o. 0.0.0.0... .048. Tsur= 24.0 
SVeETSGe wail temperature ......... Tw a= —-&. O 
initial temperature of liquid phase: Fi = 2.0 
avetage temperature of solid phase : Tsm = -3.0 


total time for Freez 


CONDUCTIVITY For salid phase at Tem 
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compare with these conductivities at Tsm= -3.0 C 


metal : ORR = 
Ice - £itCe = 
SETLES [ keer = 
paraiiel : kpar = 
empirical <« kemp 

geometrical: kaeao = 
Veinberg : RVeaR= 


energy denfity 
performance ferctor: 
storage factar 


eri. OO 
2. 39 
4.07 

104.18 

13°73 
3. 96 
$2.95 


to7 7 fo 


Ws mK 
WS eS 
W/mK 
W/m 43S 
N/m HK 
iW / meK 
W/ me 
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OOO 


percentage of phasechange 
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2.900 32.000 4.000 5.00qg 6.a00 7.000 


Lime-min. 
freezing vertical 


&.dg00 9.000 
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Fauperigernt Bo at 


metal > AL HAT MU 

BPOTOSITY > Pha = O. S02 
description «© BEADS, diameter: temm 
WO. (ares Doves oe gt, OOOO mm 
density Doo hs Gs YOR, O ka/mkte 
Epecific heat: cm = u8a0. 0 Usk ageK 


MELTING in a VERTICAL culinder 


temperature of heath oo... : Fb = BO. 9 
ambient temperstere oo... 0.0... 00.. : Tsur= 24.0 
average wall temperature o.....0.. : Twa= 21.9 
initial temperature of solid phase : Ti = 0.9 
average temperature of liquid phase: Tim = 10.9 
total tawe for megtane: tt = 6.0 min 


CONDUCTIVITY fur liguid phase at Tim = 10.9 C 


hae § Fre gat Chey feb oy EGA 
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i ho eas dil, “Ads bb tes Ke 
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compare with these conductivities at Tim= 10.7 C 


metal : Rats = PO. OO W/mkK 

watet [ Rhelo = 0.58 W/inkK 

SeTLesS ( ksepr = 1,09 W/mkiK 

baraliel ( trap = $00, 23 W/m 

empirical «© hemp = 7. AS W/m eK 
geometrical: kgeao = 1.78 W/mK 

Veinberg : Rveans 3.569 W/mkK 

energy density © Qf fee OQ. LFb6etO9 J/mtES 
nerformance Fector: pF = 8. 3i W/meK 


storage factar : 8 = OO. 16detiO JEW/mEESEK 
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metal 


poresity 
description 
Yol. /area 
density 
specific 


Experiment # 


26 


ALUMINUM 


phi 


Vea 


rhom= 


heat: cm 


BEADS: 


> O, 532 
diameter: 


2. OOOO mm 


2702.0 kg/m##2 
880.0 U/kg kK 


MELTING in a HORIZONTAL cylinder 


temperature of bath 


ambient 
average 
initial 
aveTrage 


temperature 


wall temperature 
temperature of solid phase 
temperature of liquid phase: 


total time for melting: 


CONDUCTIVITY 


keff= 


13.97 W/mK 
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compare with these conductivities at 


metal 
water 
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paraliel 
empirical 
geometrical: 
Veinberg 
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energy density 
performance factor: 
storage factor 
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For liquid phase at Tim 
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time-“min. 
melting horizontal 
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metal 
porosity 
description 
YO. fares 
density 


Specific heat: 


FREEZING itt a 


temperature o 


ambient tempereture 
average wail temperature 
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iment # <2? 


COPPER 
ehio= 
oP IRALS: 


3h C= 
Cif = 


0. 740 


). Sam X 4mm 

We = 0, S000 mm 

8933.0 kgo/mtt2 
377.2 JUSkORK 


VERTICAL cylinder 


@ BSUh ..ce cae seges - Th = -15. 


total time for Freezing: 


CONDUCTIVITY 


For solid phase at Tsm 


MEEKER KEBREEEREEEREREREER EERE 


+ ROS ine oS. 


ii W/mkK 
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initial temperature of liquid phase: Ti = 
average temperature of solid phase : Tsm = -4. 
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compare with these conductivities at Tsm= -4.8 C 


metal 

ace 

SETIES 
parallel 
empirical 
geometrical: 
Veinberg 


db 
oO 
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n 
MWe: 


sn 

kKice 
koer 
kpat 
kemp 
ggeo 
EVES 


Ho on on OR 


energy density = Ae 
performance factor: pF 
storage FfactarTr 7 8 
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O. 307etO9 J/mkkG 
6.19 W/ mith 
O. 190et1O JEW/meEgEK 


percentage of phasechange 
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Experiment # 28 


metal COPPER 

porosity phi = 9. 940 
description SPIRALS; 0. Smm X 4mm 
Vol. /area va = 9. 5000 mm 
density rhom= 8933.0 kg/mk#2 
specific heat: cm = 377.2 J/kgtK 


MELTING in a VERTICAL cylinder 


temperature of bath 
ambient temperature 
average 
initial 
average 


wall temperature 
temperature of solid phase 
temperature of liquid phase: 


total time for melting: 


CONDUCTIVITY for liquid phase at Tim 


kefF= 


3.01 W/mkK 


ee 2©« e©« * * &€ e@ @ © 
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kim= 1.49 W/m 
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compare with these conductivities at 


metal Km = 
water kh2o = 
series kser = 
parallel kpar = 
empirical kemp = 
geometrical: kgeo = 
Veinberg kvein= 


energy density. 
performance factor: 
storage factor 
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metal 
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description 
WOl. /aTres 


density 


Specific heat: 


Experiment # 


COPPER 


phi = 


SPIRALS: 


Ve 
rhoam= 
th = 


temperature of bath 
ambient 
average 
initial 
average 


temperature 
wail temperature 
temperature of solid phase 
temperature of liquid phase: 
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a? 


O. FAO 


0. Sam X 4mm 
= 0, 5000 mm 


8933.0 kg/me#2 


377, 2 J/kakk 


MELTING in a HORIZONTAL cylinder 
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total time For melting: 


CONDUCTIVITY For liquid phase at Tim 
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1.41 H/meK 
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compare with these conductivities at 


metal 
water 


SETI E5 
parallel 


empirical 
geometrical: 
Veinberg 


energy density 
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Experiment # 30 


description 
Yol. /area 


density ; 
epecific heat: 


FREEZING in a VERTICAL 


CARTRIDGE BRASS 


pha = 0. 9796 

7O% Cu, BO* Zni SPIRALS 
ve = 0.0484 mm 

Thoam= G530.0 ka/mtte 
cm = 374.6 J/kgtK 


temperature oF bath 


ambient 
average 
initial 
average 


total time for 


CONDUCTIVITY for 


temperature 


temperature of liquid phase: 
temperature cf solid phase 


cylinder 
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cclid phase at Tsm 


compare with these conductivities at Tsm= 
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energy density 
performance factor: 
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Experiment # dai 


metal : CARTRIDGE BRASS 
porosity > phi = 0.956 

description : 70% Cu, 3O% Zni SPIRALS 
¥WOl. fares : ¥e = 9.0484 mm 

density > Them= 8530.0 kg/mtte 
specific heat: cm = 374.5 J/kg tK 


MELTING in a VERTICAL culinder 


temperature oF bath ............0.. : Fb = 8O.9 

ambient temperature ...........2.. : Tsur= 22.0 

average wall temperature ......... : Twa= 26.6 

initial temperature cf solid phase : Ti = -O.7 

average temperature of liquid phase: Tim = 13.3 

total time For melting: tt = 28.0 min 

CONDUCTIVITY for liquid phase at Tim = 13.3 C 
kefF= fe. fF W/mK 

HEEHHEREKEKRKKKERRHKEREEEE 

# k Limes 1,619 W/meK + 

KRHHEKKERKERKERERREEKEERHREEE 

compare with these conductivities at Tim= 13.9 ¢ 

metal : km = 105. 95 W/mK 

water >: kheo = 0.58 W/mtK 

SeTies : kser = O. 61 W/mK 

parallel : &kpar = 3. 265 W/mtK 

empirical : kemp = 0.74 W/mK 

geometrical: kgeo = 0. 62 W/mK 

Veinberg © kvein= 0.57 W/mK 

energy density » 77% 9. 334etO9 J/mke3 


performance factor: pF = 2. 49 W/mK 
storage factor : 5 = . B382etO9 JEW/meEE4 KK 
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Experiment # 32 


metal 
porosity 


description 
vol. /area 


density 


specific heat: 


CARTRIDGE BRASS 


phi = 


704 C 
Ve = 
Thom= 


0. 756 
vu, 30% Zn 
0.0484 mm 


SPIRALS 


8530.0 kg/mteo 


€ i = 


374.6 J/kg#K 


MELTING in a HORIZONTAL cylinder 


temperature of bath 
ambient temperature 
average wall temperature 
initial temperature of solid phase 


average temperature of liquid phase: 
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total time for melting: 


CONDUCTIVITY For liquid phase at Tim 
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2. 73 W/mkK 


ee e©¢ «© * © @® oe 8 @ 


KHKKHEKKRERKKRERKERRRKEREERERKRER 


# 


kim= 


1.23 W/mekK * 
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compare with these conductivities at 


metal 
water 
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Experiment # 33 


metal : STAINLESS STEEL 

porosity : phi = 0. 959 

description : AISI 304; SPIRALS; O. imm 
Vol. /area : ve = 9.090484 mm 

density : Them 7900.0 ke/mt te 


specific heat: cm $56.8 J/kgkK 


FREEZING in a VERTICAL cylinder 


temperature of bath .............. : Jb = -15.90 
ambient temperature .............. >: Fsur= 24,90 
average wall temperature ......... >: Fwsa= -11.9 
initial temperature of liquid phase: Ti = 0.3 
average temperature of solid phase : Tsm = -6.0 
total time for Freezing: tt = 56.9 min 
CONDUCTIVITY For solid phase at Tsm = -6.0 C 
KEKHEE EEE EEE REE RR REE RRR 

* k Sin= 2. 34 W/mkK 
KERR EKKRER ERE KERERREE 

compare with these conductivities at Tsm= -6.0 C 
metal : km = 34.28 W/mK 

ice > kice = 2.20 W/mK 

S5€T2esS >: kser = 2.28 W/mK 

parallel : Kkpar = 2.569 W/mK 

empirical : kemp = 2. 37 W/mK 

geometrical: kgeo = 2. 30 W/mK 

Veinberg : kvein= 2. 38 W/mK 

energy density © @777= 9. BSIidetO9 JU/meeo 
performance factor: pF = 2. 50 W/mK 


storage factor O. Bi4etO9 JHW/meEa KK 


HOI 


a 


echande 


age of phas 


ar 
i, 


percen 


eal 


10. 


Of - 


.O0 - 


40 


00 


mG 


1} C0 


10. 


a 


g 


a .00 


STEEL 





30.00 40.00 506.00 60.00 70.00 


fime“min. 
freezing vertical 


Bo .aG 


= .00 


1 


G 


OG. 


9LT 


metal 
porosity 


description 
Vol. /area 


density 


specific heat 
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Experiment # 34 


STAINLESS STEEL 


O. imm 


phi = 0.959 

AISI 304; SPIRALS: 

ve = 9. 0484 mm 

Tham= 7900.0 ka/mtte2 
cm = $56.8 J/kgth 


MELTING in a VERTICAL cylinder 


temperature oF bath 
ambient temperature 


*s ee 8© 8&® 8 «© © » © #® @ oe #8 @ 
eo 28 @© © @© © © 8 © © #@ 8 @ © 


average wall temperature 


initial temperature cf solid phase 


average temperature Gof liquid phase: 


total time for melting: 


CONDUCTIVITY 


ke FF= 


1.97 W/meK 


HHREHHKKEREREREERERRERERRERER 
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compare with these conductivities at 


metal 
water 
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parallel 


empirical 
geometrical: 


Veinberg 


energy density 
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Experiment # 35 


metal >: STAINLESS STEEL 

POTOSI TY >: phi = 0.959 

description <: AISI 304; SPIRALS; O. imm 
YGOl. /ares : ve = 0.0484 mm 

density : Thom 7900.0 ka/m#tte2 
specific heat: cm = $56.8 JU/kqgktK 


MELTING in a HORIZONTAL cylinder 


temperature of bath .............. : Fb = GO. 
ambient temperature ......00........ >: Fsur= 24. 
average wall temperature ......... : Ju.a= 26. 
initial temperature of solid phase : Ti = -—@. 
average temperature of liquid phase: Tim = 12. 
total time for melting: tt = 23.0 min 
CONDUCTIVITY For Ligvid phase at Tim = i3a.2 C 
kefF= 2. HO W/mkK 


HREHEKKEKEKEKREKEEKKERERERREEK 


# kim= 1.08 W/meK 
KRKHEKKRRKR KHER KREKEEEEEREE 


compare with these conductivities at Tim= i322 
metal > okm = 14.28 W/mK 

water : kKheco = 0. 58 W/mkK 

Series > kser = 0. 51 W/mtkK 

parallel : kper = 1.15 W/mK 

empirical : kemp = 0.57 W/mK 

geometrical: kgeo = 0. 62 W/mK 

Veinberg : kvein= 0.55 W/mK 

energy density : Q2°7?= QQ. 3B5etO9 J/mktko 
performance factor: pF = 3. O07 W/ mith 


storage factor :  § O. 1L0Get1O JHEW/MmEE4EK 
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Appendix B. Derivations 


Appendix B includes the 


5252.3 Sud 53220248. 


As long as the fusion 
frozen fraction is 

F 

The total volume for phase 


B) in a triangular unit of 
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derivations of Eq. 5.2.2.2, 

fronts are not in contact the 
_ V(t) 
Tae 

change V (see Fig. Appendix 


all 
the tube bank is the area of the 





triangle VA times the depth 1 without the volume of the 
tubes V : The areas related to V. and V are shown in 
tube A all 
Fig. App-.B. ° 
22h ete O° tee oe a 
A 2 mie 4 
1 
‘Subs ~ 2 eae 
V = V - V 
all A tube 
The actual frozen volume is 
eee 
V(t) = > 0 t1 3 R wl 
Substituting V(t) and Vou into the definition for the 
frozen fraction 
2 2 
F = V(t) _ r -R 
A 5 
“ra. ote 9-4 2 
L ae « 
21 


and solving for r 
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ght j ee = 1/2 
r= R C[—52—-(@) - 1] F +1 ) (5.2.2.3) 
L 
The fusion fronts meet at r = 9 thus 
7 + 12 7 R2 
PB aee = 1/2 (Sie 2e2e2) 
3 2 2 
——_ L —-R 
27 


After the fusion fronts make contact with one another, the 
calculation is continued for inward freezing. The frozen 


fraction is now 


2 
2 r 
=. V 
oe a de (*) 
all V 
The total volume M ag inside the virtual cylinder of radius 
R is 
V 
7 2 
V a Rul 


and the actual frozen volume is 


- 2 2 
V = R«wt1l- r rl 
Vv Vv 


where r_ is the radius of the assumed cylindrical fusion 
Vv : 
front which contains the same volume as the actual unfrozen 


part in the heat exchanger. 


2 


rom Vol V(t) 


Substituting the last three equations into (*) and solving 


for r yields 
V 


183 


; 1/2 
R n1-(V -V(t)) 
Vv all 
r= Rk (1 - 35+ —____) 
i ¥ R-tl 
y 1/2 
Vall v(t) “all 
RY ( 2 ~ ¥ 2 ) 
R «wl all Ro«l 
Vv Vv 
o 1/2 
all — 
= RL (— [1 - Fl] ) 
Roni 


substituting V yields finally Eq. 5.2.2.4. 
a 


11 


Figure B.l Volumes in a triangular unit of staggered 
tubes. 


